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Resonant X-Ray Scattering from the Quadrupolar Ordering Phase of 

CeB 6 
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We theoretically investigate the origin of the resonant x-ray scattering (RXS) signal near the 
Ce Liu absorption edge in the quadrupolar ordering phase of CeBg, considering the intersite 
interaction between the Ts states in the initial state. The anisotropic charge distribution of the 
4/ states modulates the 5d states through the intra-atomic Coulomb interaction and thereby 
generates a large RXS superlattice intensity. The temperature and magnetic field dependence 
indicates that the induced dipolar and octupolar orders have little influence on the RXS spectra, 
in good agreement with the recent experiment. 
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Resonant x-ray scattering (RXS) has attracted much 
interest as a useful tool to investigate the orbital or- 
der, since the resonant enhancement on the superlattice 
Bragg spots for the orbiial order was observed near the 
Mn K edge in LaMnC>3.EP Although similar experiments 
have already been carried out on several materials, the 
origin of the RXS signal is still controversial. The RXS 
is described by a second order process that photons are 
virtually absorbed by exciting a core electron to unoccu- 
pied states and then emitted by recombining the excited 
electron with the core hole. Therefore, the Ap states 
are involved near the K edge for LaMnC>3 in the dipo- 
lar process. In the early stage, the Coulomb interaction 
between the 3d and Ap states was considered as the ori- 
gin J3# but subsejauejit; studies based on the band struc- 
ture calculations^' BETH' have revealed that the Coulomb 
effect is two order of magnitude smaller than the effect of 
the Jahn- Teller distortion, indicating that the Ap states 
are sensitive to electronic structure of neighboring sites. 

Recently the RXS experiment was carried out near 
the Ce Liu absorption edge in the quadrupolar ordering 
phase of CeBg.Er The lattice distortion is not observed 
because of much weaker coupling between lattice and 
well-localized 4/ states. Also, the 5d states are consid- 
ered to be more localized than the Ap states in LaMnC>3. 
Therefore, the mechanism of the RXS in CeB-6 is ex- 
pected to be different from that in LaMnC>3. The pur- 
pose of this paper is to make clear the origin and to 
discuss the implication of the RXS spectra in compar- 
ison with the experiment, by considering the intersite 
interaction between the r§ states. 

In CeBg, each Ce atom is considered to be in the f 1 - 
configuration. Within the J = 5/2 subspace, the Tg 
quartet states have lower energy than the T7 doublet un- 
der the cubic crystal field. The IVI^ excitation energy 
is ~ 540 K. With decreasing temperatures the antiferro 
quadrupole (AFQ) order appears at Tq — 3.3 K with an 



ordering wave vector Q = (1/2, 1/2, 1/2). Antiferromag- 
netic phases exist at even lower temperatures, but will 
not be considered in this paper. These phase transitions 
originate from the intersite interaction, which lifts the 
degeneracy of the T$ states. 

The intersite interaction was derived on the basis of a 
RI£KX interaction by Ohkawa.oEP Recently, Shiina et 



al.i3c3t3 extended his model by fully taking account 
of the symmetry of the interaction as well as the order 
parameters. Thereby, they solved a longstanding contro- 
versy between the neutron diffraction!!^ and NMRliS 1 in 
the context of the induced order parameters under the 
external magnetic field. Following their study, we use 
in the initial state the Hamiltonian H of the so called 
extended Ohkawa modelx3 

H = D J2 [(1 + *)M< ' M,- + t?tV + ecr, ■ a, 
(id) 
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where fi, r', ry, £ acting on the Tg states are explic- 
itly defined in ref. 13. Since the RXS spectra are not 
sensitive to the parameter values of S and e, we simply 
assume 6 — 0.2 and e = 1. This choice is krpspi to. give 
a reasonable description for the AFQ phaseOl!3'E3 As- 
suming D = 1.83 K, we get the transition temperature 
Tq = 3.3 K for the AFQ phase with no external magnetic 
field in the mean field approximation. 

Near the Ce L\\\ absorption edge, the 2p core electron 
is virtually excited to the bd states and is recombined 
with the core hole in the dipolar process. Since the 2p 
states are well localized around Ce sites, the scattering 
tensor can be approximated by a sum of the contribu- 
tions from each site of the created core hole. The corre- 
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sponding cross section is given by 
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Here G is the scattering vector, and u> is the frequency 
of photon. The j runs over Ce sites. The cross section 
becomes order N, the number of Ce sites. The dipole op- 
erators x a (jYs are defined as X\{j) = x, X2(j) = y, and 
Xs(j) = z in the coordinate frame fixed to the crystal 
axes with the origin located at the center of site j. The 
energy eigenvalues are defined as E n (j) (n = 1 ~ 4) and 
the eigenfunctions as \ipn(j)) a ^ s ^ e 3 m the mean field 
approximation. We take thermal average on these states 
with probability p n {j). State |A) represents the inter- 
mediate state with energy E A . The life-time broadening 
width r of the core hole is assumed to be 2 eV. We study 
only the situation of G = (1/2, 1/2, 1/2) in the following. 
For this G value in the conventional scattering geome- 
try^ the geometrical factors P M for incident photons 
and P' M for scattered photons are given bycS* (P a )i = 
(P^') 1 = (cos/3cosip + smip)/V2, (P (T ) 2 = (P'^'h = 
(cos/3cos-0-sin-0)/\/2, (P CT ) 3 = (P /CT ') 3 = -sin/3cos-0, 
(P ,7r ')i = [-sin0(cos/3sinV> - cos^) + cos0sin/?]/\/2, 
(P /7r ') 2 = [-sin6»(cos/3sinV> + cosVO + cos#sin/?]/V2, 
(P ,7r )3 = [sin 6 sin f3 sin -0 + cos 8 cos /3] , with 6* being the 
Bragg angle and /3 = arccosl/\/3- The azimuthal angle 
ijj is defined such that the crystal axis (1,-1,0) lies on 
the scattering plane at ip = 0. 

In the intermediate state, the 5d states of Ce_are in- 
volved. From the band calculation for LaBgjLZP it is 
known that the 5d states of La hybridize with B 2p states, 
forming energy bands with their width ~ 5 eV. The 5d 
states in CeBg are considered to be similar to those in 
LaB6, except for the relative position to the 4/ levels. 
Note that the 5d states are more localized than the 4p 
states in LaMn03, since the latter band width is as large 
as ~ 20 eV. Therefore, the dominant contribution to 
the RXS spectra comes from the local process which can 
be expressed in terms of the local Green's function con- 
structed by the 5d density of states (DOS). We roughly 
simulate it by a simple rectangular shape with the width 
5 (2) eV and the edge located 2.5 (5.5) eV above the 4/ 
levels, for the t^ g (e g ) symmetry. The hd band is slightly 
filled by electrons, whictL mediate the intersite interac- 
tion mentioned aboveBliliE 3 ! In the present calculation, 
however, we simply assume the vacant 5d band in the 
initial state, since the actual filling modifies little the 
RXS spectra. Thus, the local Green's function for the 
5d states is defined by 
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dicating either ti g or e g states. 

The intermediate state at the site of core hole are spec- 
ified by m d s d (s d : spin) of the 5d electron, m* s^ {nv = 
-3, . . . , 3) of the 4/ electron, and m p s p (m p = -1,0, 1) 
of the 2p hole. Then the resolvent in the space of these 
states is given by 
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where p 5d d (e) is the DOS mentioned above with m d in- 



where H[ nt represents the Hamiltonian in the intermedi- 
ate state. The e/ and e p represent the average energy 
of the 4/ electron and 2p hole. The interaction V con- 
tains the spin-orbit interaction for each particles, and the 
intra-atomic Coulomb interaction between them. The 
parameter values for them are evaluated in the 2p 5 4f 1 5d 1 
configuration of_the Ce 3+ atom within the Hartree-Fock 
approximationES 1 It is known that the Slater integrals 
P°(5d,4/), F°(2p,4f) and F°(2p,5d) are well screened 
while the other integrals are not much even in solids. We 
assume relatively small values, 4 eV, 12 eV, and 7 eV 
for P°(5d,4/), P°(2p,4/) and F°{2p,5d), respectively, 
while we use the atomic values multiplied by a factor 
0.8 for other quantities. Note that the anisotropic parts 
F\bd, 4/), F\bd, 4/), G 1 (5d, 4/), G 3 (5d, 4/), G 5 (5d, 4/) 
in the Coulomb interaction between the 5d and 4f elec- 
trons, not the isotropic part P°(5d, 4/), are relevant to 
giving rise to the RXS intensity. The crystal field en- 
ergy, which is important for selecting the r 8 states as 
the initial state, is negligibly small in the calculation of 
the intermediate state. The right hand side of eq. (5) 
is a matrix with dimensions 840 x 840, which we nu- 
merically invert. This resolvent is the same at all sites 
of core hole, so that the scattering amplitudes become 
different at different sublattices of the AFQ order only 
after multiplying the matrix elements of the dipole oper- 
ators between the initial and the intermediate state (see 
eq. (3)). 

Now we discuss the calculated results. We have three 
types of possible ordered phase for H = 0, in which one of 
the staggered quadrupole moments, (O xy ) (= 4(r y cr z )), 
(O yz ) (= 4(t v ct x )), and (O zx ) (= A{T y a y )), is finite. Here 
(X) indicates the thermal average of operator X. We 
simply call them as the O xy , O yz , and O zx phases. As- 
suming three domains with equal relative volumes, we 
average the RXS intensities over the domains. Figure 
1 shows the RXS spectrum as a function of photon en- 
ergy, at T = 2.7 K and ip = 0. We obtain a single 
peak structure in agreement with the experiment, u 1 
We adjusted the core hole energy e p such that the peak 
is located at the experimental position fiui = 5722 eV. 
The w-dependence is the same in the three domains. We 
roughly estimate that the peakJntensity is less than 10~ 2 
or 10~ 3 of that for LaMn0 3 MLl ) The contribution of the 
quadrupolar process is included in the curve.H Although 
it is expected to give a peak in the pre-edge region, it is 
too small to be seen. This contrasts with the RXS spec- 
tra from the AFQ phase of DyB2C2, where a pre-edge 
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peak is conspicuous .EJElP The core hole potential can- 
not become a main origin for the RXS spectra, since it 
works the same way at all sites. This point is clarified by 
comparing the spectra with those calculated by turning 
off the potential. The latter also consists of a single peak 
located at the position about 19 eV higher than the for- 
mer. In the figure, the spectrum is shifted such that the 
peak is located at the same position; the spectral shape 
and intensity are similar except for the energy shift. This 
behavior is quite different from the absorption spectra, 
which arc known to be considerably enhanced by the core 
hole potential. 
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Fig. 1. RXS spectra as a function of photon energy for G = 
(1/2, 1/2, 1/2). H = 0, T = 2.7 K, and i/i = 0. The contributions 
from three domains are averaged. The solid and dotted lines 
represent the intensities for the a — > a' channel and the a — » 
it' channel, respectively. Curves with letter "a" represent the 
spectra calculated by turning off the core hole potential. 



We obtain the scattering amplitudes A1(G,oj) in sim- 



forms, 
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in the O xy , O yz , and O zx phases, respectively. Therefore, 
the principal axes, for example in the O xy phase, are 
given by rotating the x, y axes with 7r/4 around the z 
axis. Note that the eigenstates with the eigenvalue 1 (- 
1) for the operator O xy has a larger charge distribution 
along the (1,-1,0) ((1,1,0)) direction. The principal 
axes are consistent with the anisotropic charge density 
of the 4/ states in the AFQ phase. This contrasts with 
the situation of LaMnC>3, where the neighboring oxygen 
potentials -afe-.more important through the Jahn- Teller 
distortion.™^ 

Figure 2 shows the azimuthal angle dependence of the 
intensity of the peak at hui = 5722 eV. The contributions 
from the three domains are separately shown. The de- 
pendence for the a — > a' channel is quite different from 



that for the a 
O yz and O 



channel. The curves for the O 



xy 



phases can be transformed into those for 
and O xy phases, respectively, by shifting ip 
with 27r/3. This threefold symmetry around the (1, 1, 1) 
direction perfectly matches the relation between the or- 



theOj^, O 



der parameters of the three domains. 

Under the external magnetic field, the staggered oc- 
tupole moment and the staggered dipole moment are in- 
duced in addition to the staggered quadrupole moment. 
To make clear the effects of such induced moments on the 
RXS spectra, we calculate the RXS intensity of the peak 
at Tllo = 5722 eV as a function of H . The temperature 
is set rather low, T — 1.37 K, so that the quadrupole 
order parameter is sufficiently developed for all values of 
H (the possible antiferromagnetic phase is disregarded). 
Figure 3 shows the calculated results in comparison with 
the order parameters. For H || (001), the primary or- 
der parameter is the staggered quadrupole moment of 



(O xy ). The staggered octupole moment (T x 



')) 



is induced with increasing H , but the staggered dipole 
moment is not. As shown in the left panel, the RXS in- 
tensity changes little with changing H, indicating that 
the induced octupole moment is irrelevant to the RXS 
spectra. For H | (110), the primary order parameter is 
the staggered quadrupole moment (O yz + O zx ) /v2, and 
the staggered dipole moment (J z ) (= (fycr 2 + ^r z a z )) 
is induced in addition to the induced octupole moment 
with increasing H. As shown in the right panel, the 
RXS intensity changes little with increasing H, indicat- 
ing that the induced dipole moment is also irrelevant to 
the RXS spectra. This marks the difference between the 
RXS and the neutron scattering experimentp-the latter 
of which detects the induced dipole moment 13 
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Fig. 2. Azimuthal angle dependence of the RXS intensity of the 
peak at hui = 5722 eV, for (a) the a — > a' channel and (b) the 
a — * it' channel. T = 2.7 K and H = 0. The solid, broken, and 
dotted lines represent the contributions from the Oxy, O yz , and 
O zx phases, respectively. 



Figure 4 shows the temperature dependence of the 
RXS intensity of the peak at hiu = 5722 eV, for H = 
and for several values of.ff || (1,1,-2). The inset shows 
the staggered order parameter ((O xy ) in the O xy phase) 
for H = and its square, calculated within the mean 
field approximation. The temperature dependence of the 
spectra for H = looks similar to that of the square of 
the order parameter. The dependence changes little with 
changing .values of H, in agreement with the recent ex- 
periment .eP 

In summary, we have calculated the RXS spectra near 
the Ce L\n absorption edge, considering the intersite in- 
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Fig. 3. RXS intensity of the peak at %ui = 5722 eV as a function 
of H. Left panel: H || (001); Right panel: H || (110). T = 1.37 
K and ip = 0. The solid and dotted lines represent the intensities 
for the cr — > ct' channel and the ct — > 7r' channel, respectively. 
Upper panels show the staggered quadrupole moments, (0 X y), 
(Oyz), (Ozx), the staggered dipole moments (J z ), and the stag- 
gered octupole moment (T xyz ). 
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Fig. 4. Temperature dependence of the RXS intensity of the peak 
at hui = 5722 cV, for H = and for finite H || (1, 1, -2). The 
curve for H = is not the limit of // — > but the domain 
average. Inset: the staggered quadrupole moment and its square 
for H = 0. 



teraction between the r§ states in the initial state. We 
have demonstrated that the 4/ charge distribution mod- 
ulates the 5d states through the intra-atomic Coulomb 
interaction and thereby generates a large RXS super- 
lattice intensity. This mechanism is different from that 
in LaMnOa, where the oxygen potentials at neighbor- 
ing sites are much important. We have evaluated the 
temperature and magnetic field dependence of the RXS 
intensity, in good agreement with the recent experiment. 
The magnetic field dependence indicates that the in- 
duced dipolar and octupolar orders have little influence 
on the RXS intensity. Thus, it can be concluded that 
the RXS spectrum in CeBg is a direct reflection of the 
AFQ order, in contrast with the neutron diffraction and 
NMR experiments. We have used a simplified model for 
the 5d bands. A change in the 5d DOS shape and the 
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filling actually modifies the spectral intensity, but does 
not alter the qualitative features discussed in this paper, 
such as a single-peak shape, the temperature and the 
magnetic field dependences. A quantitative study using 
a more realistic model is left for future. 
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